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ABSTRACT 

The "algorithm driven by the density estimate for the identification of clusters" (DED- 
ICA, Pisani 1993, 1996) is applied to the A3558 cluster complex in order to find 
substructures. This complex, located at the center of the Shapley Concentration su- 
percluster, is a chain formed by the AGO clusters A3556, A3558 and A3562 and the 
two poor clusters SC 1327-312 and SC 1329-313. We find a large number of clumps, 
indicating that strong dynamical processes are active. In particular, it is necessary to 
use a fully three-dimensional sample (i.e. using the galaxy velocity as third coordi- 
nate) in order to recover also the clumps superimposed along the line of sight. Even 
if a great number of detected substructures were already found in a previous analysis 
(Bardelli et al. 1998), this method is more efficient and faster when compared with 
the use of a wide battery of tests and permits the direct estimate of the detection 
significance. Almost all subclusters previously detected by the wavelet analyses found 
in the literature are recognized by DEDICA. 

On the basis of the substructure analysis, we also briefly discuss the origin of 
the A3558 complex by comparing two hypotheses: 1) the structure is a cluster-cluster 
collision seen just after the first core-core encounter; 2) this complex is the result of a 
series of incoherent group-group and cluster-group mergings, focused in that region by 
the presence of the surrounding supercluster. We studied the fraction of blue galaxies 
in the detected substructures and found that the bluest groups reside between A3562 
and A3558, i.e. in the expected position in the scenario of the cluster-cluster collision. 

Key words: galaxies- clusters- substructures- individuals: A3558- A3562- A3556- 
SC 1329-313- SC 1327-312- SC 1329-314 



1 INTRODUCTION 

Clusters of galaxies are thought to form by accretion of sub- 
units in a hierarchical bottom up scenario. This view arises 
naturally in various theories of cosmic structure formation, 
as f.i. the Cold Dark Matter dominated scenarios. Numerical 
simulations on scales of cosmological relevance (e.g. Cen & 
Ostriker 1994, Colberg et al. 1997) revealed that mergings 
happen along preferential directions, called density caus- 
tics, which define matter flows, at whose intersection rich 
clusters are formed. Detailed high resolution simulations of 
such collisions were performed by McGlynn & Fabian (1984) 
and Roettiger, Burns & Loken (1993), who studied cluster- 
cluster and cluster-group merging, respectively. In the sim- 
ulations, the two cores have two or three encounters before 
dissipating all the kinetic energy: in particular, the less mas- 
sive cluster emerges as a partially dispersed galaxy popula- 
tion along the down merger side of the dominant cluster. 
Members of the merging cluster, trapped nearby the dom- 



inant core, create an apparent dumpiness on either side of 
the center. The kinetic energy is dissipated mostly by the 
particles in external regions of the clusters, which expand 
creating a common envelope, while the cores do not change 
dramatically. 

From the observational point of view, the detection of a 
large fraction of clusters which present substructures (see f.i. 
Kriessler & Beers 1997 and references therein) revealed that 
these systems are still dynamically young and show evidence 
of mergings. The best studied examples are A2256 (Briel, 
Henry & Bohringer 1992), where a small group is detected 
in the X-ray band nearby the cluster center, and Coma, 
where a number of substructures are revealed (Biviano et al. 
1996). In particular. Burns et al. (1994), on the basis of the 
detection of a large scale X-ray emitting filament connecting 
Coma with the NGC4839 group, suggest that the group has 
probably already passed the cluster center. 

The most spectacular phenomena are expected to be 
seen when the merging involves repeated collisions or even 
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coalescence of comparable units: during the intermediate 
stages, the clusters form structures and complexes on scales 
of few megaparsecs. 

Rich superclusters are the ideal environment for the de- 
tection of cluster mergings, because the peculiar velocities 
induced by the enhanced local density of the large scale 
structure favour the cluster-cluster and cluster-group col- 
lisions, in the same way as the caustics seen in the simu- 
lations. The most remarkable examples of cluster merging 
seen at an intermediate stage is found in the central re- 
gion of the Shapley Concentration, the richest supercluster 
of clusters found within 300 h~^ Mpc (Zucca et al. 1993; 
hereafter h=Ho/100). On the basis of the two-dimensional 
distribution of galaxies of the COSMOS/UKSTJ catalogue 
(Yentis et al. 1992), it is possible to find several complexes of 
interacting clusters, which suggest that the entire region is 
dynamically active. Therefore, this supercluster represents a 
unique laboratory where it is possible to follow cluster merg- 
ings and to test related astrophysical consequences, as the 
formation of radio halos, relicts and wide angle tail radio 
sources and the presence of post-starburst (E-l-A) galaxies. 

In this paper we perform a substructure analysis of the 
A3558 cluster complex using a non-parametric, scale inde- 
pendent algorithm, based on the adaptive kernel density 
estimator. This method could be considered alternative or 
complementary to the wavelet decomposition analysis, ex- 
tensively used in literature. In particular, a wavelet analysis 
of this complex was done by Girardi et al. (1997) and, lim- 
ited to the dominant cluster A3558, by Dantas et al. (1997): 
these analyses offer an opportunity to test the two methods 
on a real, rather complex situation. In Section 2 we briefly 
describe the method, in Section 3 the used samples are pre- 
sented and in Section 4 a brief discussion of the subclumps 
is reported; finally in Section 5 we present our conclusions. 



kernel basically because it is differentiable in all its domain. 
This is a crucial feature in order to identify the density peaks 
according to Eq. ^. Even if the most efficient (i.e. with min- 
imum variance) kernel form is the Epanechnikov function, 
Silverman (1986) showed that difference in the use of the 
two forms is negligible. Indeed, what is crucial for cluster 
detection is the width and not the shape of the kernels, pro- 
vided that they satisfy the normalization and convergence 
condition (see section I of Pisani 1993) . 

The estimate of the kernel sizes is done through an it- 
erative procedure, starting from a large value for the kernel 
sizes and progressively reducing them until the minimum of 
the "cross-validation function" is reached. The various steps 
are the following (for details see Pisani 1993, 1996): 
1) set the first guess (n — 0) of the size as ctq — 4at where 



at = 0.96Af' 



-l/(d+4)l 



\ 



(3) 



where sf is the standard deviation of the coordinate (Sil- 
verman 1986). In this step, a large starting ao is selected: 
note that for a symmetric gaussian distribution with disper- 
sion ao and same number of points of our samples, we have 
(Jo ~ 1 — I.Sctg; 



2) set the following iteration (n ) as a„ 



i/2; 



3) take as "pilot function" fp{x) the density kernel estimate 
fk (x) with sizes ai fixed at an Vi (see eq. 27 and 28 of Pisani 
1993); 

4) set a'i = \ia„, where Ai is a "local bandwidth" defined as 



where 



(4) 



2 THE METHOD 

In order to identify the substructures in our sample and to 
assess their significance, we followed the method known as 
"algorithm driven by the density estimate for the identifica- 
tion of clusters" (hereafter DEDICA) developed by Pisani 
(1993, 1996). This is a non-parametric, scale independent 
method, whose output gives a list of groups with the re- 
lated statistical significance, their members, and a list of 
isolated objects. Here we remind briefly the main steps of 
the method. 

First, we estimated the density field f{x) (where x is the 
n-dimensional position vector) using the "adaptive kernel 
probability estimator" for which 



1 ^ 

fk{x) = —^K{x^,a^;x) 



(1) 



where fk{x) is the kernel estimator of f{x), Xi is the position 
of the i*'' galaxy, A'^ is the number of objects in the sample 
and ai is the size of the gaussian kernel defined as 



K{xi, ai; x) 



1 



(27ra2) 



2W/2 



exp 



2af 



(2) 



where d is the number of dimensions of the sample. 

We stress that we choose the Gaussian form for the 



1 

log 9 = j^^^og{fp{xi)) 



(5) 



Now, the sizes of the kernels a'i are depending somehow from 
the local density through the factor Ai: the higher is the 
density, the smaller are the a'i\ 

5) evaluate the cross-validation function CV[fk\{an) of the 
n*'^ estimate of the density fk (eq. 4 and 5 of Pisani 1996). 
The CV[fk\{an) estimates the quadratic deviation between 
fk{x) and the "true" parent density function f{x) (see eq. 
33 of Pisani 1993). 

The steps from (2) to (5) are iterated until the minimum of 
CV[fk\{a7i) is reached: the corresponding kernel widths a'i 
give the optimal density estimate fk{x) as: 



1 

fk{x) ^ K{xi, ai; x) 



(6) 



Peaks in the density distribution function are assumed 
to define the presence of clusters (in this section we use 
the word cluster in the statistical sense): the local maxima 
of fk{x) are found by solving for each starting i^^ galaxy 
position {i — 1,N) the iterative equation 



Xm-\-l,i — Xm.i ~\~ 



EN 
i = l 



fk{Xm. 



^fk{Xm.i) 
fk ( ^ m , i ) 



(7) 
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where xi^i is the original position of the i"' galaxy. Eq. (Q) 
defines a path from the original position to a limiting point 
along the maximum gradient of the function fk (x) All start- 
ing points which reach the same limiting position xnm for 
m — > oo define a cluster. From the analysis of the distribu- 
tion of the distances between the various XHm in our samples, 
we decided to consider coincident all the limiting points that 
differ by less than 0.1 arcmin. 

The significance of the cluster is estimated by calculat- 
ing first the likelihood function of the sample 



L = l[hix,) = l[J2^^^^ 



(8) 



= 1 M = 



where Ff_i(xi) is the contribution to the density function of 
the z^*** cluster C^, calculated as 



(9) 



where the sum is done on the galaxies belonging to the /x*'' 
cluster and u is the number of clusters found in the sample. 
For /i = we refer to the set of isolated galaxies, i.e. those 
that did not reach a xum shared with other object. 

Then the likelihood of the /x*'' cluster (I/^j) is calcu- 
lated by replacing in Eq. the sizes a[ of the kernels of 
the substructure objects with abck, i.e. the kernel size of the 

background. The ratio = ~21n -p- is distributed as a chi 
square variable with one degree of freedom and is related 
to the probability that the density distribution is better de- 
scribed by the presence of the /x*'' cluster rather than by a 
uniform background. 

In an ideal case, it would be possible to determine atck 
by using the isolated galaxies, because they are likely to 
be part of the background. However, also local fiuctuations 
of the density field of a cluster can be identified as local 
maxima and hence this fact causes a misidentification of 
cluster members as isolated galaxies. 

A method to identify really isolated galaxies is to check 
the relative contribution to the total density of the kernel of 
the i*'' galaxy in the position Xi (Pisani 1996), given by 



P{i G Co) 



1 K(xi,ai;xi) 



N 



fk{x) 



(10) 



This equation estimates the probability that a "nominally" 
isolated galaxy is really part of the background. We found 
that for our samples P{i £ Co) < 0.15, hence we can con- 
clude that there are no really isolated galaxies. Note that 
the estimate of P{i G Co) is independent from the cluster 
identification step (Eq. 

In order to determine a value for the background, we 
plotted the histogram of the kernel sizes a'^, which is a sym- 
metric distribution with a tail at higher values. We choose 
as atck the value at the beginning of the tail (see Fig. |^ in 
Sect. 3.1). 

The group significances calculated from the values of 
are strictly correct only if the background is well determined, 
that is not our case. For this reason we regard the only as 
a ranking parameter for the existence of the groups and we 
consider only subclusters in the tail at higher values of this 
parameter (see Fig. ^ in Sect. 3.1). This choice corresponds 
approximately to a "formal" significance of 99.9%. We forced 



the background at various reasonable densities in order to 
check the robustness of our results, finding that the general 
clustering pattern of our samples is not strongly dependent 
on the poorly known global background density. 

As a test of the efficiency of this method, we applied 
DEDICA to simulated Gaussian distributions, where no sub- 
structures are expected. We performed various simulations 
of 2-D and 3-D Gaussians, with a number of points between 
300 and 700 (which are the approximate numbers of galax- 
ies in our various samples), finding that the success rate 
(i.e. absence of significant substructures) of DEDICA is al- 
ways between 80% and 90% (depending on the number of 
objects). In the cases in which DEDICA finds significant 
substructures, the numbers of these "spurious" groups is of 
the order of 2-3. 

This method differs from the wavelet decomposition be- 
cause it is non-parametric and scale-independent. Indeed, 
the wavelet method is a deconvolution procedure, which as- 
sumes an exploring function, which is fixed throughout the 
whole sample volume. On the contrary, DEDICA adopts a 
kernel function, which is locally density dependent and in 
this sense is not a deconvolution method (Silverman 1986). 
The wavelet analysis is able to detect structures depending 
on the scale of the exploring function, while DEDICA does 
not requires any assumption on the scale of clusters. 



3 THE A3558 CLUSTER COMPLEX 

In a series of papers (Bardelli et al. 1994, 1996, 1998, here- 
after Paper 1, 2 and 3; Venturi et al. 1997) we performed 
a detailed multiwavelength study of the A3558 cluster com- 
plex, a remarkable association formed by three AGO (Abell, 
Corwin & Olowin 1989) clusters A3556, A3558 and A3562, 
already noticed by Shapley (1930). This structure is approx- 
imately at the geometrical center of the Shapley Concentra- 
tion and can be considered the core of this supercluster. By 
the use of multifibre spectroscopy, we obtained 714 redshifts 
of galaxies in this region, confirming that the complex is 
a single connected structure, elongated perpendicularly to 
the line of sight. In particular, the number of measured red- 
shifts of galaxies belonging to A3558 is 307, and thus this 
is one of the best sampled galaxy clusters in the literature. 
Moreover, two smaller groups, dubbed SC 1327-312 and SC 
1329-313, were found both in optical (Paper 1) and X-ray 
band (Breen et al. 1994; Paper 2). In particular, in Paper 2 
we detected a bridge of X-ray emitting hot gas connecting 
A3558 and SC 1327-313 (and possibly extending to A3562, 
see figure 7 of Ettori et al. 1997) which resembles the struc- 
ture between Coma and NGC4839 described by Burns et al. 
(1994). The presence of this feature strongly supports the 
suggestion by Metcalfe et al. (1994) that A3562 already en- 
countered A3558: the bridge could be the tail of gas left by 
A3562 after the first passage through A3558 (for details see 
Burns et al. 1994). 

The clusters of this complex attracted the attention of 
various investigators, who applied several statistical tests in 
order to detect substructures. Bird (1993), applying both 
2-D and 3-D tests to galaxies in A3558, found three sub- 
clumps: the first coincides with the cluster center and cor- 
responds to the main component, the second lies eastward 
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of the cluster center and the third is in the south-east quad- 
rant. 

In Paper 3, applying the test developed by Dressier & 
Shectman (1988), we detected a group at ~ 10 arcmin from 
the position of the second clump of Bird (1993). From the 
analysis of the isodensity contours and the velocity distribu- 
tion histograms, there was evidence of another density excess 
at ~ 15 arcmin southward of the A3558 center. Moreover, 
A3558, SC1329-313 and A3556 appear to have bimodal ve- 
locity distributions. 

Girardi et al. (1997), using a modified version of the bi- 
dimensional wavelet decomposition which takes into account 
also the dynamical information, individuated in A3558 the 
same subclumps found by Bird (1993) and by us (Paper 3), 
and an additional substructure westward of the center of 
A3562. 

Dantas et al. (1997), using the bi-dimensional wavelet 
decomposition in A3558, detected a bimodality in the core 
of the cluster. The authors claim that this result resolves 
both the cD offset problem (in fact the dominant galaxy 
was suspected not to be at rest with respect to the cluster 
velocity centroid, Gebhardt & Beers 1991), and the "/3 prob- 
lem" (Paper 2), by decreasing the central velocity dispersion 
from ~ 1000 km/s to ~ 800 km/s 

Analysing a ROSAT PSPC image of A3558, we noted 
(Paper 2) that the diffuse emission of the hot gas is better 
described by two components, one centered on the geomet- 
rical center of the cluster and the other on the cD galaxy. 
Slezak et al. (1994) applied the wavelet analysis to this X- 
ray map and found evidence of a substructure nearby the 
cluster core. 



3.1 The hi— dimensional sample 

From the COSMOS/UKST galaxy catalogue (Yentis et al. 
1992), we extracted a rectangular region of 3°. 2 x 1°.4, 
corresponding to 13''22'"06" < a(2000) < 13''37™15' and 
-32°22'40" < (5(2000) < -30°59'30". This region is part of 
the UKSTJ plate 444 and contains the A3558 complex. We 
restricted our analysis to the 2241 galaxies with magnitudes 
brighter than bj = 19.5. The isodensity contours of galaxies 
of this region are presented in Fig. |^. 

The adopted ai,ck is 6.4 arcmin, which corresponds to 
0.05 gal arcmin"'^, and a threshold of 12 (Fig. □). 

The DEDICA algorithm found 113 groups, 16 of which 
are considered significant. In Table |l| the 16 significant clus- 
ters sorted by right ascension are listed. Column (1) is the 
identification number while columns (2) and (3) give the a 
and S coordinates of the group center. We identified these 
points as the common xum position of the members and 
they do not necessarily coincide with the geometrical cen- 
ters. Columns (4) and (5) report the number of the clus- 
ter members and of the available velocities in the group. 
Columns (6) and (7) are the dynamical parameters com- 
puted using the biweight estimators of position and scale 
(see Paper 1 for details on the method). Finally, Column 
(8) reports the main component to which the substructures 
belong, i.e. the A3556, A3562, A3558 Abell clusters, and 
the poor clusters SC 1317-312 and SC1329-313 (denoted as 
"Poor"). 

The velocity sample (see below) covers an area which 
is ~ 17% of the bi-dimensional sample and therefore not all 




13» 35"> 13i> 30" 13" 25" 

Right Ascension 

Figure 1. Galaxy isodensity contours of the region of the A3558 
cluster complex: the objects are binned in 2 arcmin X 2 arcmin 
pixels and then smoothed with a Gaussian with a FWHM of 6 
arcmin. Dashed lines correspond to 1 Abell radius circles around 
cluster centers. The positions of the two poor clusters SC 1327- 
312 and SC 1329-313 are shown. 



subclusters have measured velocities: for this reason in two 
cases the dynamical parameters are left blank in Table ^ 



3.2 The three— dimensional sample 

In this case we used the velocity sample described in Paper 
3. From the 714 galaxies with measured redshift we selected 
the 600 objects in the velocity range [10000 - 20000] km/s. 
The velocity sampling of the A3558 complex is not constant: 
it varies from ~ 20% in the region of A3562 to ~ 48% in 
the core of A3558. This fact could cause a bias that erases 
undersampled substructures or leads to false detections in 
regions more sampled than the average. 

In order to check this point, we performed the bi- 
dimensional analysis on this subsample. Among the 14 bi- 
dimensional groups present in the area of the velocity sam- 
ple, 13 are correctly detected and B2240 (corresponding to 
T599) has been lost because of a significance lower than the 
limit. The tight correspondence between bi-dimensional and 
these "three-dimensional" groups indicates that this bias is 
not severe. 

Given the fact that the velocities are not related to po- 
sitions because of the dominance of peculiar motions over 
the Hubble fiow, a problematic point is how to handle the 
velocity coordinate. Our formalism uses symmetric three- 
dimensional gaussian kernels, with = ay = cr„ and there- 
fore it is necessary to scale the velocity interval in order to 
have a numerical range comparable to the other two vari- 
ables. The typical a of the spatial coordinates is 2-3 arcmin, 
while the typical cr„ (determined on the one-dimensional red- 
shift distribution) is ~ 200 km/s. We have therefore cho- 
sen to compress the velocity by a factor 100; we have also 
checked that the method gives essentially the same results 
if this reduction factor is in the range 50 — 500. 

A more precise approach to this problem requires the 
modelling of the peculiar velocity field. This implies the de- 
tailed knowledge of the mass distribution in the A3558 com- 
plex, which is presently not available. In fact there is no 
general agreement even for the mass of the dominant, best 
studied, cluster of the complex (A3558). The mass obtained 
from optical data ranges from 3.4 x 10^* h"^ M© (Dantas 
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Table 1. Groups in the bi-dimcnsional sample 



# 



a (2000) 5 (2000) # mem. # vel. 



<v > 



notes 



B109 


13 


22 


31.4 


-31 


18 


22 


34 




B194 


13 


22 


57.7 


-31 


37 


46 


67 


14 


B353 


13 


24 


15.4 


-31 


40 


40 


46 


22 


B694 


13 


26 


58.3 


-31 


21 


29 


57 


20 


B830 


13 


27 


07.2 


-31 


51 


53 


72 


26 


B906 


13 


27 


52.0 


-31 


29 


23 


81 


52 


BIOOO 


13 


28 


04.0 


-31 


44 


41 


55 


24 


B1014 


13 


28 


05.8 


-31 


05 


24 


24 




B1056 


13 


28 


07.4 


-31 


33 


05 


53 


29 


B1199 


13 


29 


00.5 


-31 


49 


23 


51 


18 


B1311 


13 


29 


23.7 


-31 


21 


34 


57 


23 


B1443 


13 


30 


28.0 


-31 


36 


41 


100 


37 


B1573 


13 


31 


08.2 


-31 


46 


58 


111 


49 


B1952 


13 


33 


28.1 


-31 


41 


16 


179 


35 


B2103 


13 


34 


50.8 


-31 


40 


14 


37 


8 


B2240 


13 


35 


38.3 


-31 


55 


40 


76 


18 



14360: 
14492" 



-117 
-145 
-131 
-164 



140381224 

1472811^,^ 

138131^36 

140641112 
14728^279 



15029 

14652 
14059 



+228 
-375 



-182 
-149 
-116 
-257 



1 4447+230 
-305 

147511^^^ 
14343lt« 



520 
575 



115 

68 
-t-86 

44 



674ll« 
lOOllii^ 

991+178 
-122 

10811^35 



842_ 
1051 



-t-150 
56 

-f 160 
-53 



1113 



1462 



103 
112 
-1-396 
-107 
+445 
100 



A5556 
A3556 
A3556 

A3558 
A3558 
A3558 
A3558 
A3558 
A3558 
A3558 
A3558 

Poor 
Poor 

A3562 
A3562 
A3562 



et al. 1997) to 6 x lO" h'^ Mq (Biviano et al. 1993), while 
the X-ray data give mass in the range 1.5 — 6.0 x lO^'' h"^ 
Mq (Paper 2, Ettori et al. 1997). 

The adopted value of the background is abck ~ 7.5 ar- 
cmin, corresponding to 0.003 gal arcmin"'^ (km/s/100)~i 
and a significance threshold of = 20. The significant 3- 
D clusters are 21, out of a total of 36, and are reported in 
Table ^ Column (1) is the subcluster identification number. 
Columns (2), (3) and (4) report the coordinates (a, S and 
Vo) of the group center, while the number of members is 
reported in Column (5). The estimated dynamical parame- 
ters are reported in Columns (6) and (7) and in Column (8) 
the general position is given. Note that in some cases, the 
velocity distributions are highly asymmetric and therefore 
there is a significant difference between Vo (the position of 
the density peak) and the mean velocity estimated with the 
biweight method. 

In Figs. 1^, ^, 1^ and ^ the positions of the substructure 
members are shown, superimposed on the isodensity con- 
tours of galaxies obtained with a smoothing of 2 arcmin. 
Note that in these figures the isodensity contours are ob- 
tained directly by binning the data in one arcmin cells and 
not from the ft'- however the two methods give similar con- 
tours. In Figs. §, I, I and |l^, the groups members are shown 
projected on the X- Velocity plane. Solid lines connect the 
position of each galaxy to the common limiting position xum 
of the group. 

4 RESULTS AND DISCUSSION 
4.1 A3556 

The bi-dimensional analysis found 2 significant groups in 
the inner part of this cluster, being the substructure B109 
at about one Abell radius from the center. The distribution 



of galaxies is not smooth as can be seen also in Fig. ^ where 
the isodensity contours show three peaks. However, the dy- 
namical analysis of the substructures reveals that mean ve- 
locities and velocity dispersions are consistent with the val- 
ues of < u >= I4357I76 km/s and a = 643^43 km/s given 
in Paper 3 for the whole cluster. 

The three-dimensional analysis shows four subclumps: 
two of them (T49 and T96) have dynamical parameters con- 
sistent each other at about one sigma and correspond to 
two nearby (~ 12 arcmin) density maxima separated on the 
plane of the sky. 

The group r75 appears superimposed along the line of 
sight with T96 (they correspond to B353), but at higher 
velocity: the presence of these two peaks in the redshift 
distribution was already detected in Paper 3, with consis- 
tent dynamical parameters. Moreover, as noted by Venturi 
et al. (1997, 1998), r75 contains the head-tail radiosource 
J1324 — 3138: the head-tail appearance of this source could 
indicate an interaction with a gaseous medium during the 
infall of the group toward the main clump of A3556. 

The group ri7 seems to be at an intermediate mean 
velocity and separated in the plane of sky from the other 
clumps. The westernmost galaxy of this group (well outside 
the density peak) corresponds to the Wide Angle Tail radio 
source studied by Venturi et al. (1997). In the bi-dimensional 
analysis this object is associated to B109, but also in this 
case it is at the edge of the substructure. This fact is anoma- 
lous for Wide Angle Tail sources, found normally at the cen- 
ter of clusters or groups (Owen 1996). 

In conclusion A3556, which appears at a first look as 
a relatively relaxed cluster (see figures 1 and 2 of Paper 3), 
revealed after a more detailed analysis a relevant dumpiness 
in its core. This is an indication that the dynamical processes 
have not had enough time to destroy the substructures and 
to smooth the galaxy distribution. 



© 0000 RAS, MNRAS 000, 000-000 



6 S. Bardelli et al. 

Table 2. Three-dimensional sample 
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4.2 A3558 

The bi-dimensional analysis revealed that the inner core 
(r < 0.75 h~^ Mpc ) of this cluster could be separated in 
two groups: _B906, which contains the dominant galaxy, and 
-B1056. This spatial bimodality was already found by Dantas 
et al. (1997) using the bi-dimensional wavelet analysis: the 
positions of their clumps A and B differ only by ~ 1 arcmin 
from our groups and there is a good consistency in the dy- 
namical parameters. However, the three-dimensional analy- 
sis shows that the situation is more complex, with four sub- 
clusters (TSST, T281, T300 and T338) almost aligned along 
the line of sight. In particular, the bi-dimensional group 
S906 is broken into two components, at < v >= 15456 
km/s (r281) and < v >= 13597 km/s (r337): note that the 
mean velocity of T281 is consistent with that of clump A' 
of Dantas et al. (1997). 

The main three-dimensional structure is the group 
T337. The group T194 has similar dynamical parameters 
and could therefore be considered an elongation of T337. 
Merging these two clumps together in order to have a larger 
sample, we find < v >= 13640lj24 km/s and cr = 665lgg 
km/s. We are consistent with the results found by Bird 
(1994), Girardi et al. (1997) and Dantas et al. (1997) that the 
velocity dispersion of the main component of A3558 is of the 
order of ~ 700 km/s. This velocity dispersion is more con- 
sistent with the lower X-ray temperature found by ROSAT 
(3.25 kcV Paper 2), rather than the higher temperature de- 
tected by ASCA (5.5 keV, Markevitch & VikhUnin 1997). 
The velocity histogram of T337 -|- T194 is asymmetric with 
a tail toward low velocities: the peak of the distribution is 



located at ~ 14100 km/s, a value consistent with the cD 
velocity. 

The bi-dimcnsional clump BIOOO has a positional coin- 
cidence in the plane of the sky with T301. These structures 
are located on a density excess, southward of the cluster 
center, already noted in Paper 3 and detected by the use 
of the West et al. (1988) symmetry test. This density ex- 
cess was revealed also by the presence of two peaks in the 
velocity distribution in the third radial shell of Paper 3. 
The dynamical parameters of the lower velocity peak are 
well consistent with those of r301, while the higher veloc- 
ity peak is probably associated to galaxies of group T413 
which also fail in this shell. The latter clump has been de- 
tected by Bird (1994), Girardi et al. (1997) and in Paper 3 
using the Dressier & Shectman (1988) A test. As noted in 
Paper 3, this substructure is in the region of A3558 which 
show a significant enhancement of temperature (Markevitch 
& Vikhlinin 1997). 

Another relevant condensation is the group T260 (or 
B830), in the south-west part of A3558, whose isodensity 
contours appear elongated. 



4.3 SC1327-312 and SC1329 -313 

An interesting feature of the A3558 complex is the presence 
of a number of clumps between A3558 and A3562, which 
form a series connecting the two clusters, starting with T413 
and T359 just on the east of A3558. In the X-ray band there 
are two relevant diffuse emissions dubbed SC 1327-312 and 
SC 1329-314 (Breen et al. 1994; Paper 2) already noted also 



© 0000 RAS, MNRAS 000, 000-000 



A substructure analysis of the ASS 5 8 cluster complex 7 







hi 


i 








a) - 


- 








I 


\ 




- 




1 















2 4 6 8 10 

o"V [arcmin] 



o 



6 

3 




Figure 2. a) Histogram of the kernel sizes cr^ for the bi— 
dimensional sample. The vertical dashed line corresponds to the 
adopted crtcfc = 6.4 arcmin. b) Histogram of the values for 
the detected groups in the bi— dimensional sample. The vertical 
dashed line corresponds to the adopted threshold of significance 
(X^ = 12). 

in the optical band (Paper 1). The SC 1327-312 group prob- 
ably corresponds to T478, even if there is an offset of ~ 5 
arcmin between the two positions. In Paper 3 we noted a bi- 
modality in the velocity distribution of SC 1329-313: these 
two groups are recovered as T520 and T496. 

4.4 A3562 

Prom the analysis of the three-dimensional sample it is pos- 
sible to find out three components. The main cluster com- 




20 40 60 80 100 

X (arcmin) 



Figure 3. Groups in the A3556 cluster region. The positions of 
sub-cluster galaxies found in the three-dimensional sample are 
overplotted to the 2-D isodensity contours smoothed with a 2 ar- 
cmin FWHM Gaussian. Symbols used to label the different groups 
are explained in the figure. 

ponent appears to be T561, with a velocity dispersion of 

~ 900 km/s, consistent with that found in Paper 3. Note 
the presence of other two clumps (T599 and T598) proba- 
bly connected with the main cluster. 



5 CONCLUSIONS AND DISCUSSION 

The DEDICA algorithm is a non-parametric, scale inde- 
pendent method for substructure analysis. We used this al- 
gorithm on both a bi-dimensional and a three-dimensional 
(i.e. using the galaxy velocity as third coordinate) sample 
in the region of the A3558 complex. The bi-dimensional 
analysis finds a large number of clumps in correspondence 
of evident density peaks, but it appears insufficient to cor- 
rectly characterize all the substructures: in fact the three- 
dimensional analysis finds several cases in which different 
subclumps aligned along the line of sight are seen as single 
bi-dimensional peaks. 

On the other hand, the not uniform and incomplete 
sampling in velocity may in principle lead to miss some 
structures in the three-dimensional analysis. However, this 
incompleteness does not seem to introduce significant arti- 
facts, because there is a good correspondence between the 
three-dimensional subclustcr centers and most of the bi- 
dimensional density excesses. 

A number of these subclusters were already detected in 
Paper 3 by the simultaneous use of various techniques, such 
as the classical shape estimators, the West et al. (1988) and 
the Dressier & Shcctnian (1988) tests and the analysis of 
the velocity histograms with the KMM algorithm (Ashman, 
Bird & Zepf 1994). The efficiency of these methods depends 
on the properties of the substructures: for example, the West 
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Figure 4. Projection on the X- Velocity plane of groups in the 
A3556 cluster region. Solid lines connect the position of each 
galaxy to the common limiting position of the group. 

et al. test recovers groups as departures from the mirror sym- 
metry of the main cluster, while the Dressier & Shectman 
method is insensitive to subclusters which are superimposed 
along the line of sight or with similar mean velocity and 
velocity dispersion. This fact led Pinkney et al. (1996) to 
recommend for the substructure analysis of clusters the use 
of a wide battery of statistical tests. Under this aspect, the 
DEDICA algorithm is more efficient and faster and permits 
a direct estimate of the significance of the detected clump, 
without relying on long computer simulations. 

A direct comparison with the results obtained with the 
wavelet decomposition analysis is not straightforward be- 
cause of the different starting samples. The subclusters de- 
tected by Girardi et al. (1997) are a subset of our sample, 
revealing that the two methods are similar as discussed by 
Fadda et al. (1998). The difference in the dynamical param- 
eters are probably due to the fact that the method of Gi- 
rardi et al. is not fully three-dimensional. Our bi-dimensional 
analysis of the center of A3558 gives the same results of 
Dantas et al. (1997), i.e. the inner 0.5 Abell radius region is 
formed by two clumps, one of which (_B906) approximatively 
centered (within ~ 1 arcmin) on the X-ray excess of the cD 
galaxy. On the contrary, our three-dimensional analysis re- 
veals a more complicated situation. 

It is possible to consider two main hypotheses for the 
origin of the A3558 cluster complex; one is that we are ob- 
serving a cluster-cluster collision, just after the first core- 
core encounter, while the other considers repeated, incoher- 
ent group-group and group-cluster mergings focused at the 
geometrical center of a large scale density excess (the Shap- 
ley Concentration). 

In the first scenario, we could speculate that a clus- 
ter collided onto A3558 and its remnants are visible as the 
overdensity regions of SC 1327-312, SC1329-313 and A3562: 
the connection of hot gas between A3558 and SC1327-312 
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X (arcmin) 




-40 -20 20 40 
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Figure 5. The same as Fig. for groups in the A3558 cluster 
region. For clearness the group members have been split in two 
panels 

found in Paper 2 could be a trace of such a collision. In 
this framework, A3556 would be formed by the members of 
the intervening cluster remained in the backside part of the 
merging direction. Therefore, galaxies outside A3558 would 
belong to the destroyed cluster and to the outer regions of 
the more massive target cluster and form the dumpiness 
expected in the simulations. 

Caldwell & Rose (1997) suggested that the infall of a 
group onto a cluster may trigger starburst phenomena on 
its galaxies. Indeed they found a number of objects showing 
recent star formation in the filament connecting the Coma 
cluster with NGC4839, which was suspected by Burns et al. 
(1994) to have already passed the core of Coma. In order to 
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Figure 6. The same as Fig. |4| for groups in the A3558 cluster 
region. 



Figure 8. The same as Fig. ^ for groups in the Poor cluster 
region. 




check if the cluster interaction modifies the galaxy spectral 
emission, we cross-correlated our velocity sample with the 
photometric sample of Metcalfe et al. (1994), which contains 
B magnitudes and the {U — B), {B — R) and {R — I) colours 
to _B = 20.8. The use of the velocity sample, limited in 
the [10000 — 20000] km/s range, eliminates the ambiguity 
of possible colour efl^ects due to objects at various redshifts. 
We consider "blue galaxies" those with ([/ — B) < 0.3: this 
choice should allow to exclude the "red sequence" in the 



colour-magnitude planes typical of the presence of clusters 
and to have at the same time a large number of objects (157) 
in the blue sample. 

Dividing the galaxies into two groups, i.e. those that are 
in significant groups and those that are outside them, we 
found that the fractions of blue objects are 25% and 45%, 
respectively. The reddest group is T337 with a fraction of 
5% of blue galaxies and this supports the idea that the group 
could be considered as a preexisting, already relaxed cluster 
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Figure 10. The same as Fig. ^ for groups in the A3562 cluster 
region. 



(that is A3558). The bluest groups are r520 and r441 with 
47% and 45% of blue galaxies and are in the expected posi- 
tion in the scenario of a cluster-cluster collision, i.e. between 
A3562 and A3558. 

In the second scenario, instead of a single merging phe- 
nomenon, there would have been a series of minor merg- 
ings and collisions, confined in a relatively small region by 
the deep potential well of the supercluster. In this case the 
different "colours" of the groups would be related only to 
the different initial galaxy morphological composition of the 
merging components. 

The determination of the density excess in galaxies 
in the entire supercluster (Bardelli et al., in preparation), 
through the analysis of a inter-cluster galaxy redshift sur- 
vey (Bardelli et al. 1997) , and of the implied peculiar veloc- 
ities, will help in discriminating among these two scenarios: 
in fact, peculiar velocities of the order of 1000 km/s at the 
center of the Shapley Concentration would point toward the 
hypothesis of a cluster-cluster collision, while significantly 
lower velocities would be consistent with the second sce- 
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